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Abstract 

Dynamical correlations of the spin and the energy density are investigated in 
the critical region of the random transverse-field Ising chain by numerically 
exact calculations in large finite systems {L < 128). The spin-spin autocorre- 
lation function is found to decay proportional to (log t)-2a.m g^j-^^j ^jQg ^^-2x'^ 
in the bulk and on the surface, respectively, with Xm and x^ the bulk and 
surface magnetization exponents, respectively. On the other hand the critical 
energy-energy autocorrelation functions have a power law decay, which are 
characterized by novel critical exponents r]e ~ 2.2 in the bulk and ??| ~ 2.5 
at the surface, respectively. The numerical results are compared with the 

predictions of a scaling theory. 
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The asymptotic behavior of the time-dependent correlation functions for interacting 
many-body systems turned out a very difficult subject of theoretical research. Exact re- 
sults in this field are scarce, one can mention the one-dimensional spin 1/2 XY-model |1| 
and the Ising chain in a transverse- field 0. Both models can be mapped onto a system of 
non-interacting fermions, where the equal-position correlation functions are calculated by 
the Pfaffian method utilizing the theory of Toplitz determinants. 

In this Letter we consider - at the first time - the critical dynamical correlations of an 
interacting quantum system in the presence of quenched (i.e. time-dependent) disorder. 
It has recently become clear that quenched disorder has rather different effects on phase 
transitions in quantum systems than on those thermally driven phase transitions. For 
example, in the Griffiths phase, which is situated at the disordered side of the critical point, 
the susceptibility has an essential singularity in classical systems, whereas in a quantum 
system the corresponding singularity is stronger, it is in a power law form. 

Here we consider the prototype of random quantum systems the one-dimensional random 
transverse-field Ising model defined by the Hamiltonian: 

H = -Y: Jicrfaf^, - E hicrf , (1) 
I I 

where the cxf , af are Pauli matrices at site / and the Ji exchange couplings and the hi 
transverse-fields are random variables with distributions tt{J) and p{h), respectively. The 
Hamiltonian in (|l|) is closely related to the transfer matrix of a classical two-dimensional 
layered Ising model, which was first introduced and studied by McCoy and Wu [^. 

The static critical behavior of the random transverse-field Ising model in (|ID has been 
studied analytically and numerically by several authors The system possesses a crit- 

ical point at 5 = [In J]av — [ln/i]av = 0, and has a spontaneous ferromagnetic order if the 
average couplings are stronger than the average fields. (We use the bracket [. . .]av to denote 
disorder averages.) The critical properties of the model, which are known through exact 
and conjectured results to a large extent, are in many respects different from that of pure 
systems. One important difference, that in the random system - due to a broad distribution 
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of various physical quantities - the typical and average quantities are usually different and 
the rare events dominate the critical properties. For instance the static average spin-spin 
correlation function is expected to behave as 

GTir) = [{afaf^r)] = ^ exp(-r/0 , (2) 

L J av rpLd^rn 

where (. . .) means the (zero-temperature) expectation value. For the random transverse- 
field Ising model the average correlation length ^ ~ b^^ diverges with the true exponent 
z/ = 2 and the decay exponent Xm = 1 — uj/2 ^ .191 is expressed in terms of the golden 
mean cu = (1 + ^/5)/2. The decay of the average end-to-end distance critical correlations 
involves the surface magnetization exponent = 1/2. On the other hand the typical 
correlation length diverges with z/typ = 1 and the typical critical correlations are of a stretched 
exponential form: — logG'5^p(r) ~ r^^^. In contrast the critical energy-density correlation 
function (r) = (o"fcrf_^^) is a self- averaging quantity and at the critical point it behaves 
as — log G'^(r) ~ r^/^, like its typical value. 

In this Letter we consider the time-dependent correlation functions 

GTir, t) = [{afit)af^,.)U and G^r, t) = [{afit)aU,.)U (3) 

at the critical point, both in the bulk and at the surface of the system. In a quantum 
system statics and dynamics are inherently related and the time evaluation is given via the 
Heisenberg picture by crf{t) = exp{tH)af exp{—tH). For simplicity here we confine ourselves 
to the autocorrelations, i.e. r = 0, dynamical two-site correlations will be discussed elsewhere 

i- 

To start our study we present a scaling framework for the quantum critical dynamics 
of the model (|l]). Consider the general time and position dependent correlation function 
(erf (t)crf_|_^), which can be written as 

(^rw«^r+.) = E(oi^r»(^i^/+rio) exp[-t(K - ^o)] . (4) 

n 

Here |n) denotes the n-th excited state of H in eq. (^ with energy En. Before performing the 
disorder average we note that this correlation function is not self averaging at the critical 



point. To see its scaling behavior at the critical point we present the following simple 
argument. The random samples can be divided into two groups. In the typical samples (i.e. 
which appear with probability one) the critical correlations decay faster than any power law. 
On the other hand a vanishing fraction of the samples (the so called rare events) is ordered 
at the critical point and the correlation function measured on these samples is of order 0(1). 
The disorder average of the correlation function is then determined by the rare events and 
the corresponding scaling behavior is governed by the scaling properties of the probability 
distribution of these rare realizations. 

For example the probability P{1), which measures the occurrence of samples with a finite 
local magnetization m{l) = 0(1) at site / (take for instance fixed boundary conditions, or 
consider an off-diagonal matrix element in the case of free b.c. , see p[), scales as the average 
critical magnetization P{l/b) = b^^"^P{l), when lengths are rescaled by a factor 6 > 1. For 
equal time correlations in the rare realizations the local magnetization is of order 0(1) 
at both spatial coordinates. The corresponding joint probability distribution P2{l,l + r) 
factorizes for large spatial separations lim^-^oo -^2(^5 1 + r) = P{l)P{l + r), since the disorder 
is uncorrelated. Consequently the spatial correlations follow the scaling rule: 

^^(r, t = 0) = r'"'"G™(r/6, t = 0) , (5) 

whereas for end-to-end distance correlations we have the surface magnetization scaling di- 
mension x^. Now taking r = b we recover the known critical decay as given in eq. (|^). 

For critical time-dependent spin-spin autocorrelations, however, the scaling behavior is 
different from that in eq. (^. This is due to the fact that the disorder is strictly correlated 
along the time axis and the probability for the occurence of a rare sample with m{l) = 0(1) 
at different times is simply P2{{l,t), {1,0)) ~ P{1)- Thus the scaling behavior of the critical 
magnetization autocorrelation function satisfies the scaling rule: 

G"^(r = 0,lnt) = r^'"G"(r = 0,\nt/b^/^) , (6) 

where we have made use of the relation between the relevant time tr and length ^ scales 
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~ In tr Note that the usual scahng combination is t/V , however, the critical 

dynamical exponent 2; is cxd here. Taking now the length scale as 6 = (Int)^, we obtain 

G™(r = 0,t)~(lnt)-2-"^ (7) 



For the surface autocorrelation function the scaling relation in eq. and consequently 

the decay exponent involves the surface magnetization exponent x^. 

For energy density autocorrelations the typical realizations govern the scaling properties 
at the critical point. The relevant quantity is now the matrix-element [I (Ol'^z^l^) P]av 

on the r.h.s. of eq. which scales in an exponential form: log[| (0|(Tf |?2)p]av = 

6^-*^/^ log[|(0|(jfyj,|n)p]av [§]• Consequently the critical energy density autocorrelations sat- 
isfy the scaling relation: 

logG'(r = 0,lnt) = fe^^/^ logG'"(r = 0,lnt/6^/2) , (8) 

and with h = (Int)^ one obtains a power law dependence of G"^(r = 0,t) with novel, non- 
trivial exponents: 

G^{r = 0,t) r-. r"^ . (9) 

In the actual calculations we transformed the model in eq. (0) into a free fermion model 
T0| , where the correlation functions are expressed by averages of fermion operators, which 



are then calculated by Wick's theorem and by the Pfaffian method |TT|]. We use free bound- 
ary conditions, in which case the most convenient representation is given in [|1^, which 
necessitates only the diagonalization of an 2L x 2L matrix. From the corresponding eigen- 
values and eigenvectors one obtains the elements of the Pfaffian, which is then evaluated 
by calculating the determinant of the corresponding antisymmetric matrix. Details of the 
calculations will be presented elsewhere 0. 

The critical properties of the random quantum spin chains are expected to be independent 
of the details of the distributions of the couplings and the fields. In this Letter we consider 
the binary distribution 7r(J) = jS{J — A) + — A~^) and h = ho, and the uniform 
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distribution 7r(J) = G(l - J)e(J) and p{h) = ho^e{ho - h)e{h). In both cases the critical 
point is at /io = 1- AH numerical data which we present below are averaged over 50000 
samples. 

First we study the critical spin-spin autocorrelation function for imaginary times t = —it 
in the bulk (i.e. at the site i = L/2) and at the surface (i.e. at site i = 1). As shown in Fig. 
la the finite lattice results fall to the same curve for logr < ^/L and the critical temporal 
decay takes place on a logarithmic scale G^^2i'^) ~ (logr)"^^™ in agreement with the scaling 
prediction (^. For surface correlations the numerical calculation is less demanding and one 
can go up to finite systems of size L = 128. As can be seen in Fig. lb in this case the 
logarithmic decay depends on the surface magnetization exponent: G^{t) ~ (logr)~^^™. 

The autocorrelation functions in real time generally have an oscillatory character. In 
the random system the average over different oscillating functions results in a complicated 
looking behaviour, as we demonstrate it for the surface autocorrelation function in Fig. 2a. 
Its Fourier transform, however, has a nice scaling character. We actually consider 

1 fOO POO 2 

= TT / / drGTit + zr) = -\{u\a^^\0)\' , (10) 

ZVr J-oo J-oo UJ 

where {uj\ is a state with an excitation energy -Eexc. — Eq = u. For small frequencies u we 
expect the finite size scaling form of xTi^) be given by 

xr (^, L) ~ xMu)/L'/') (11) 

with the scaling combination log(ct;)/L^/^ replacing log(t)/L^/^ from (|^). In Fig. 2b we show 
a corresponding scaling plot that yields a good data collapse. 

Next we turn to analyze the energy density autocorrelation function at the critical point. 
As seen on Fig. 3a the energy density autocorrelation function is described by a power law 
dependence in imaginary time r as G'i^2i'^) ~ agreement with the scaling prediction 

(|) and (^). The decay exponent r/e — 2.2 is universal, i.e. it does not depend on the type of 
the randomness. A similar power law decay is found for the surface energy autocorrelations 
in Fig. 3b, with a surface critical exponent rj^ ~ 2.5. These novel critical exponents complete 
our knowledge about the critical behavior of the random transverse-field Ising spin chain. 
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To summarize we have studied dynamical correlations at the critical point of the random 
transverse-field Ising spin chain. We showed that the magnetization autocorrelation function 
has anomalous logarithmic decay, whereas the energy-density autocorrelations decay as a 
power law with novel critical exponents. There are still many interesting aspects of the 
dynamical behavior of random quantum systems. Here we mention the dynamical properties 
in the Griffiths phase, the temperature dependent autocorrelations and the dynamical two- 
site correlations. The study of these and other related problems are in progress 
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FIGURES 
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FIG. la: Bulk spin-spin autocorrelation function G2/2(''') — [('^L/2(^)'^L/2)]av in imagi- 
nary time for various system sizes (and the uniform distribution). Note that we have chosen 
L to be odd, so that L/2 denotes the central spin. In this plot with [G^y2 (''")] "^^^^"^ on linear 
scale versus r on a logarithmic scale the infinite system size limit is expected to lay on a 
straight line as indicated. 
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FIG. lb: Same as a) for the surface spin-spin autocorrelation function G^{t) — 
[{o'i{T)ai)]a,y in imaginary time. 
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FIG. 2a: (Top) Surface spin-spin autocorrelation function G^{t) in real time for the 
binary distribution with A = 4. The data for L=64 and those shown for L=32 are exactly 
identical, although both data sets have different disorder realization. The expected l/log(t) 
behavior for the envelope indicated by the broken line is only a guide to the eye. 
FIG. 2b: (Bottom) Scaling plot of the Fourier transformed surface spin-spin autocorrelation 
function xTi^) (HO) f*^^ binary distribution and A = 4. 
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FIG. 3a: (Top): Bulk energy-energy autocorrelation function G|,/2(''') ~ 
[{^L/2i'^)^L/2)]a-v imaginary time for various system sizes (and the binary distribution, 
A = 4) in a log-log plot. The straight line has slope —2.2, which yields our estimate for the 
exponent rje. FIG. 3b: (Bottom) Same as a) for the surface energy-energy autocorrelation 
function Gi(r) = [(c'i(''')ci)]av m imaginary time. The straight hne has slope —2.5, which 
yields our estimate for the exponent rj^. 
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